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Abstract 

Kinematic equations for the motion of slowly propagating, weakly curved fronts 
in bistable media are derived. The equations generalize earlier derivations where 
algebraic relations between the normal front velocity and its curvature are as- 
sumed. Such relations do not capture the dynamics near nonequilibrium Ising-Bloch 
(NIB) bifurcations, where transitions between counterpropagating Bloch fronts may 
spontaneously occur. The kinematic equations consist of coupled integro-differential 
equations for the front curvature and the front velocity, the order parameter associ- 
ated with the NIB bifurcation. They capture the NIB bifurcation, the instabilities 
of Ising and Bloch fronts to transverse perturbations, the core structure of a spiral 
wave, and the dynamic process of spiral wave nucleation. 



1 Introduction 



Interfaces separating different equilibrium or nonequilibrium states appear in a variety of 
contexts including crystal growth, domain walls in magnetic and hydrodynamic systems, 
and reaction-diffusion fronts [1] . The global patterns that appear in these systems depend to 
a large extent on the possible occurrence of interfacial instabilities. A transverse instability 
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of the interface, for example, may lead to fingering and the formation of labyrinthine pat- 
terns [2-8]. Another instability with dramatic effects on pattern formation is the nonequi- 
librium Ising-Bloch (NIB) bifurcation [9-13]. The bifurcation, which takes a single stable 
(Ising) front to a pair of counterpropagating stable (Bloch) fronts, has been found in chemi- 
cal reactions [14-16] and in hquid crystals [17-19]. Far below the NIB bifurcation stationary 
patterns or uniform states prevail. Far beyond it, a regime of ordered traveling patterns, 
including spiral waves, exists. In the vicinity of the bifurcation disordered spatio-temporal 
patterns, involving repeated events of spiral-wave nucleation appear. This behavior, which 
we call "spiral turbulence" , has been attributed to spontaneous transitions between counter- 
propagating Bloch fronts induced by curvature variations, front interactions, and interactions 
with boundaries [6,20]. 

A common theoretical approach to studying pattern formation in interfacial systems is based 
on a geometric equation for the interface curvature (see Eqn. (2) below) [21-25]. Once the 
dependence of the normal velocity of the interface on its curvature is known the shape of 
the interface can be determined at any given time. For reaction-diffusion fronts this depen- 
dence may become particularly simple: Away from a NIB bifurcation, a linear relation is an 
excellent approximation [24,26,27]. The curvature equation, however, does not capture pos- 
sible transitions between counterpropagating fronts near a NIB bifurcation because the front 
velocity becomes an independent slow degree of freedom and can no longer algebraically be 
related to curvature [20,28,29]. 

In this paper we consider bistable media that exhibit NIB bifurcations and derive kinematic 
front equations which generalize the geometric curvature equation. The new kinematic equa- 
tions capture transitions between counterpropagating fronts, and spontaneous spiral-wave 
nucleation, a process which plays a crucial role in the onset of spiral turbulence. The equa- 
tions are: 

• An equation for the order parameter, Co, associated with the NIB bifurcation: 

^ = (ae - a)C, -pd + ^K + ^o + ^-^J KC^ds' . (1) 



• A geometric equation for the front curvature, k: 

|._(.^ + ^)C.-|/.C„d.'. (2) 



• An equation relating the normal front velocity C„, the curvature k, and the order param- 
eter, Co'. 

Cn^Co-DK. (3) 
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In these equations s is the front arclcngth, and the critical parameter value, designates the 
NIB bifurcation point. Note that Eqn. (3) cannot be regarded as a linear relation between 
the normal velocity of the front and its curvature since Cq is not a constant but a dynamical 
variable coupled to curvature through Eqn. (1). In fact, Eqns. (1) and (3) can be recast into 
a single integro-differential equation for the normal velocity (using Eqn. (2)) 

- ^ / KCr^ds' , (4) 




dCr, 

dt 



d_ 
ds 



which replaces the algebraic Cn — k relation used in earlier derivations. An algebraic Cn — n 
relation can be recovered from Eqns. (1) and (3) assuming the order parameter Co follows 
adiabatically slow curvature variations. This issue is further discussed at the end of Section 3. 

The order parameter equation (1) yields the NIB bifurcation for planar (uncurved) fronts. 
For a symmetric bistable system (70 = 0), the Ising front, Cq = 0, is stable for a > a,.- At 
a = ac the Ising front becomes unstable and a pair of Bloch fronts appears, Cq = Cq = 
(ttc — Oi)/ (3. For a nonsymmetric system (70 7^ 0) this pitchfork bifurcation unfolds into 
a saddle node bifurcation in the usual way. 

A brief account of the results to be reported here has appeared in Ref. [30]. We present in 
Section 3 a detailed derivation of the kinematic equations for a particular reaction-diffusion 
model introduced in Section 2. In Section 4 we study the kinematic equations. We analyze 
the stability of planar fronts to transverse perturbations and present numerical solutions 
describing steadily rotating spiral waves and spiral-wave nucleation induced by a transverse 
instability. We conclude in Section 5 with a discussion of our results. 



2 The reaction diffusion model 



We consider the FitzHugh-Nagumo model with a diffusing inhibitor, 

— = e~'^{u -u^ -v) + 5~^V'^u , 
dt 

— =u — aiv — ao + V V , (5) 
at 

where u and v, the activator and the inhibitor, are real scalar fields and is the Laplacian 
operator in two dimensions. The parameter Oi is chosen so that (5) describes a bistable 
medium having two stable uniform states: an up state {u^,v^) and a down state 
Ising and Bloch front solutions connect the two uniform states {u±, v±) as the spatial coordi- 
nate normal to the front goes from —00 to +00. The parameter space of interest is spanned 
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Fig. 1. The NIB front bifurcation and planar front transverse instability boundaries in the e — S 
parameter plane. The thick line is the NIB bifurcation, 5_F(e) = ij'^/e. The dashed lines are the 
boundaries for the transverse instability of Ising, Si{e), and Bloch, Sb{(-), fronts. When 6 > 6i 
{S > Sb) planar Ising (Bloch) fronts are unstable to transverse perturbations. The thin lines are 
approximations to the transverse instability boundaries obtained from a linear stability analysis of 
the kinematic equations. Parameters: ai = 4.0, oq = 0. 

by 6,5 and Oq, or alternatively by 77 = y/eS, jJ. — e/5, and Oq. Note the parity symmetry 
(w, v) — > (— —v) of (5) for oo = 0. 

The NIB bifurcation line for ao = is shown in Fig. 1. For ^ ^ 1 it is given by 5 = 5f{^) = 
Tjl/t, 01 Tj = Tjc, where rjc = ^^^3 and = ai + 1/2 [13]. The single stationary Ising front 
that exists for rj > rjc loses stability to a pair of counterpropagating Bloch fronts at 77 = rjc- 
Beyond the bifurcation {rj < rjc) a Bloch front pertaining to an up state invading a down 
state coexists with another Bloch front pertaining to a down state invading an up state. Also 
shown in Fig. 1 are the transverse instability boundaries (for = 0), 6 = 6j{e) = e/r]'^ and 
5 = 5b (e) = r]c/-\/e, for Ising and Bloch fronts respectively. Above these lines, 5 > Sj^b, planar 
fronts are unstable to transverse perturbations [6,31]. All three hnes meet at a codimension 
3 point PS: e — rj^, 5 — 1, Qq — 0. 
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3 Deriving the kinematic equations 



The derivation of the kinematic equations is described here in three steps: A transformation 
to a frame moving with the curved front (Section 3.1), derivation of the normal velocity 
equation (3) using a singular perturbation approach (Section 3.2), and derivation of the order 
parameter equation (1) (Section 3.3). In deriving the equations we assume that /i = e/5 ^ 1 
and that curvature is small, k <^ 1. Additional assumptions needed in deriving the order 
parameter equation are described in Section 3.3. 

3. 1 The moving frame 

We transform to an orthogonal coordinate system (r, s) that moves with the front, where r 
is a coordinate normal to the front and s is the arclength. We denote the position vector 
of the front by X(s,t) = (X, F), and let it coincide with the m = contour line. The unit 
vectors tangent and normal to the front are given by 



where 9{s^ t) is the angle that s makes with the x axis. A point x = (x, y) in the laboratory 
frame can be expressed as 



This gives the following relation between the laboratory coordinates (x, y) and the coordi- 
nates (s, r) in the moving frame: 



s = cos + sin Oy , 



r — — sin ^x + cos 6y , 



x = X(s,i) +rf . 



X — X{s, t) — r 



dY 
ds 



y^Y{s,t)+r 



97' 



(6) 



where we used the fact that s — dX-fds. 



With this coordinate change, partial spatial derivatives transform according to 




d _ dY d 
dx ds dr 



dX d 
ds ds ' 



dy ds dr 



d dX d 




ds ds ' 



(7) 



where 



G = (l + rn) 
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and the front curvature, is given by 

de dY d'^X dX d'^Y 



K — — 



ds ds ds^ ds ds^ 



The time derivative transforms according to 



d dr d ds d . . 



where [23,24] 



ds 
di 



^ J nCnds' , (9) 



dr 

^ = -C„. (10) 
Using (7) we find for the Laplacian operator frame [32] 

ar^ or os os os'' 



The reaction-diffusion system (5) in the moving frame is 



r]^/jIDu — u — — V + /iV'^u , 

Vv — u — Qiv — ao + V'^v , (12) 

where V and are given by (8) and (11), respectively, and we recaU that fi = e/5 and 
77 = \fe5. 



3.2 The normal velocity equation 



We study Eqns. (12) assuming <^ 1. Note that the hmit — > can be taken safely without 
departing from the immediate neighborhood of the front bifurcation ai rj = r]^- We will use 
this fact in deriving the order parameter equation. Consider the narrow front region where u 
changes on a spatial scale of order y^/I while v changes on a scale of order unity Stretching 
the normal coordinate r according to 2; = r/^fji, Eqns. (12) become 
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--,(Vv + u-a.v-a„-G'^-G^^J^O. (13) 

Expanding m and v as 

= t^O + VJiVi + + ... , 
we find at order unity the stationary front solution 

1^0 = — tanli(2;/\/2), t;o = . 

At order y^/I we find tlie equations 

^'dt dz '"^ dz ' dz'^ 



dzduo duo d\i 
Cui = vi+r)——-KG—, ^ = 0, (14) 



wliere 



Solvability of (14) yields 



rj2 

/:^— + l-3ul. (15) 



dr 3 , ^ 



where Vf = v{0,s,t) + 0(e^) is the approximately constant value of the inhibitor v in the 
narrow [(9(y^)] front core region. The first term on the right-hand-side of (16) is identified 
with the order parameter for the NIB bifurcation: Cq = ~:^^'^f- Since the normal velocity 

is Cn = —%i Eq. (16) yields the normal velocity relation (3) with D = 



3.3 The order parameter equation 



Away from the narrow front region u and v change on the same spatial scales. Letting /i ^ 
in Eqs. (12) wc obtain u — — f = 0. The relevant solutions are u = u^{v) ^ 1 — v/2 for 
r < and u = u^{v) ^ —1 — v/2 for r > (assuming oi is sufficiently large) [13]. We thus 
obtain the following free boundary problem for v: 
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7=^f-^ + Pi + P2, 




r < 0, 



r > 0, 



(17) 



where M ^ - ^ + 




K— ao + G' 




ds dv 
dt ds ' 



ds^ 



and the square brackets denote jumps of the quantities inside the brackets across the front 
at r = 0. 

We confine ourselves to nearly symmetric systems (|ao| <S 1) and to a parameter regime that 
includes the immediate vicinity of P3 and extends into the Ising regime near or below the 
transverse instability boundary 6 = 5/(e). This allows solving the free boundary problem (17) 
by expanding propagating curved front solutions as power series in c around the stationary 
planar Ising front [28,33], where c <^ 1 is the speed of a planar Bloch front solution. We 
assume weak dependence of k and v on s and achieve this by introducing the slow length 
scale S = cs and assuming X = X(S', t). This assumption dictates k = c^kq where kq ~ 
We also introduce a slow time scale T — (?t to describe deviations from steady front motion. 
We write 



00 



v(r, t,T) = V 



(0)(^) + ^cV")(r,5,t,T) , 



(18) 



n=l 



where 



V 



V 




r < 0, 
r > 0. 



Expanding 



V = Vc- c\ + c^V2 + ••• , 
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anticipating a pitchfork bifurcation, and using these expansions in (17) produces the set of 
equations 



-p("), n = l,2,3, 



(19) 



where 



.(2) 



and 



3 



•,(3) 



3 r (1) dv^ (2) 

(1) 



3771 a) 3 



+ y(r, 5,r) + aoo-(i-r^)«o 



+ 1; 



(2) 



+ V 



(3) 



5r ' 



V{r, S, T) 



dT 



G 



ds^ ^ &r ds 



(20) 
(21) 
(22) 



(23) 



In (22) we assumed = c^floo where aoo ~ C^(l)) ^iid recall that Kq — k/c^. Notice that 
1^ ~ C'(l), and P2 contributes only at orders higher than c^. 

Equations (19) should be supplemented by appropriate asymptotic conditions as r — > ±00. 
Since lim^^qioo v^^\r) — ±g~^ and Qq ~ 0{c^), the asymptotic conditions in (17) are satisfied 
by demanding 

lim v^") = n = 1, 2 , (24) 



and 



lim ^ -ooo/g^ ■ 



(25) 



We are interested in solutions of Eqns (19) at long times where they become independent of 
the fast time scale t: dv'^'^ydt — as t ^ 00. For n — 1, the stationary solution of (19) with 
the asymptotic condition (24) is [28] 

v^'\r,S,T)=vlZ,F{r), 
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where 



F{r) = (1 - gr)e«'' , r < , 
F(r) = (1 + gr)e-*'' , r > . 

Notice that v^^'' decays to zero as |r| — oo on a scale of order ~ C'(l). Since k <^ 1 we 
approximated G in Eqn. (23) hy G = (1 + r/t)"^ fa 1. 

For n = 2 we obtain the stationary solution [28] 



;(2)(r,5,T) 



,(2) 



^ (1) ^^3, 



F(r), 



(26) 



For n — 3, the stationary solution of (19) with (25) is 



^(3)(r,5,r)=-^ + 



r < 0, 
r > , 



(27) 
(28) 



where 



/1^_„3 (1) (2) , 



4g ^ ' ' 2^ 1'"=° 77c 1''=° 1''=° 2g2 ° 



-6^ ^1^ 

Application of the (no) jump condition [v^^^J^^q = leads to 



4 2(1-5-1) a^-yW dSdv^^^ 



3g 



«o + 



(29) 



where v'^^^ is evaluated at r = 0. Using the expansion (18), the integral (9), and transforming 
back to the fast variables s, t, Eqn. (29) becomes 



dvf _ V2{r)c - 77) 2(1 - S-^) d^Vf dvj 



dt 



ml 4772 ^ 3 



3g 
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where Vf{s,t) — v{0,s,t). Equation (30) coincides with the order parameter equation (1) 
once we make the following identifications: Cq — ~:^Vf, a — ccc = /? = 1/6, 
7 = ac{l - S~'^), and 70 = 2acqao. 

Algebraic Cn — k relations can be obtained from Eqns. (1) and (3) for smooth weakly curved 
fronts assuming Cq follows adiabatically slow curvature variations: 

Cn^Co-DK, (31) 



where Cq solves 

(a;e-Q;)Co-/3Co^ + 7« + 7o = 0. (32) 



Such an assumption is vahd away from the NIB bifurcation, but the condition Co ~ c <^ 1 
used in deriving Eqn. (1) no longer holds. To test how Eqns. (31) and (32) perform away 
from the NIB bifurcation we compared them with C„ — k relations obtained from the implicit 
equation 

C^ + Dk^ + (33) 



derived in [6]. Eqn. (33) is vahd at any distance from the NIB bifurcation. Figs. 2 show 
graphs of Eqns. (31) and (32) (solid curves) and of solutions of Eqn. (33) (dashed curves) 
close to and away from the NIB bifurcation. The agreement between the two approaches 
remains very good even where c > 1 (Fig. 2b). Thus the adiabatic elimination of Co away 
from the bifurcation reproduces the algebraic Cn — n relations to a very good approximation. 



4 Numerical solutions of the kinematic equations 



We study two types of solutions to the kinematic equations: steadily rotating spiral waves 
(Section 4.1), and the nucleation of a spiral- wave pair by a transverse instability (Section 
4.2). 
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Fig. 2. The Cn — n relations (31)-(32) derived from the kinematic equations (sohd curve) and 
solutions of the implicit relation (33) (dashed line), (a) Near the NIB bifurcation, e = 0.0115, both 
relations give the same result, (b) Farther from the bifurcation, e = 0.0105, the agreement is still 
good even though the kinematic equations are derived for c <C 1. Parameters: ai = 4, ao = —0.0001, 
6 = 1.063. 

4-1 Spiral waves 

Consider a "front" solution connecting the planar Bloch front, Co = Cq, k = 0, at s = —oo 



with the planar Bloch front, Co = Cq , k = 0, at s — +oo, where Cq = ±^J {ac — a)/(3, 
and we have assumed a symmetric model, ao = or 70 = 0. Fig. 3a shows such a solution 
obtained by numerically integrating (l)-(3). As demonstrated in Fig. 3b this front solution 
of the kinematic equations (l)-(3) represents a spiral-wave solution of the FitzHugh-Nagumo 
model (5). Far away from the spiral core the leading front of the spiral approaches a planar 
Bloch front pertaining to an up state invading a down state (C Cq as s ^ —00). The 
trailing front approaches a planar Bloch front pertaining to a down state invading an up state 
(C Cq as s ^ +00). The spiral core is naturally captured as the interface separating 
these two Bloch fronts. Cores of spiral waves in bistable and excitable media have also been 
studied in Refs. [35-40] using steady state free boundary formulations with linear C„ — ac 
relations. 

Fig. 4a shows a similar front solution but for an asymmetric system, 7^ or 70 7^ 0. 
Fig. 4b shows the corresponding spiral wave. The spiral tip, defined as the point of zero 
curvature, is no longer stationary, but rotates along a circle. A word of caution is needed 
here, however. The kinematic equations do not take into account interactions between front 
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Fig. 3. A front solution to the kinematic equations (l)-(3). (a) The order parameter Cq and the 
curvature n along the arclength s. (b) In the x — y plane the front solution corresponds to a rotating 
spiral wave. The shaded (light) region corresponds to an up (down) state. Parameters: oi = 4.0, 
ao = 0, e = 0.01234, 5 = 1.0. 

segments. Including such interactions would require changing the asymptotic conditions, 
v(=Foo, s, t) = v± in Eqns. (17). As a result, away from the tip (outside the frame in Fig. 4b), 
the trailing front of the spiral meets the leading front and the up state domain disappears. 



4.2 Spiral nucleation induced by a transverse instability 

Earlier numerical solutions of the FitzHugh-Nagumo model (5) revealed that an instabihty 
of a planar front solution to transverse perturbations near a NIB bifurcation can induce 
spontaneous nucleation of spiral waves followed by domain splitting [31]. The kinematic 
equations (l)-(3) capture the transverse instability of the Ising front and, to linear order 
around the codimension 3 point, P3, also the transverse instability boundary of the Bloch 
fronts. Explicit expressions for the transverse instability thresholds for a symmetric system 
(70 = 0) can readily be obtained. Let 

Co = Cq+ Co exp(o-t + iQs) + c.c. 
K, = + R exp(at + iQs) + c.c. 
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Fig. 4. A front solution to the kinematic equations for a nonsymmctric system (70 7^ 0). (a) The 
order parameter Cq and the curvature k along the arclength s. (b) In the x — y plane the front 
solution corresponds to a rotating spiral wave. The shaded (light) region corresponds to an up 
(down) state. Parameters: a\ = 4.0, ao = —0.0001, e = 0.0115, 5 = 1.0. 



where (C°, k") = (0,0) for the Ising front and (C°, k^) = {a^ - a)/ 13,0) for the Bloch 
fronts. Inserting these forms in (l)-(3) gives the following transverse instability lines, lin- 
earized around 5 — 1: 



Ising: e — rjl5, Bloch: e — r]l{3 — 25) . 



These lines are displayed in Fig. 1 (thin lines). Fig. 5 shows typical growth rates of transverse 
perturbations of wavenumber Q for Ising fronts (solid line) and Bloch fronts (dashed line) for 
the symmetric case. Note that the first wavenumber to grow as the transverse instability lines 
are traversed is Q = 0, consistent with our assumption of small curvature in the vicinities of 
(or below) these lines. 

To test whether spiral wave nucleation induced by a transverse instability is captured by 
the kinematic equations, we numerically computed the time evolution of a planar front near 
the NIB bifurcation and beyond the transverse instability boundary. Figs. 6a-d show four 
snapshots of this time evolution. The initial front pertains to an up state invading a down 
state (Co > 0). The transverse instability causes a small dent on the front to grow (Fig. 6b). 
The negative curvature then triggers the nucleation of a region along the arclength where 
the propagation direction is reversed (Cq < 0) (Fig. 6c). The pair of fronts in the kinematic 
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Fig. 5. Growth rates, cj, of transverse perturbations of wavenumber Q to uniform solutions of the 
kinematic equations (l)-(3). Sohd hne: Bloch fronts, e = 0.011, 5 = 1.08. Dashed Une: Ising fronts, 
e = 0.012, 5 = 1.2. In both cases ai = 4.0, oq = 0. 



equations that bound this region correspond to a pair of counter-rotating spiral waves in 
the FitzHugh-Nagumo equations (Fig. 6d). Fig. 6d demonstrates the equivalence of spiral 
pair nucleation in the bistable medium to "droplet" nucleation in the kinematic equations. 
The Cn — K relation for the parameter values of Figs. 6 is shown in Fig. 2a. Although the 
relation does not capture the nucleation dynamics it does provide a heuristic explanation of 
the nucleation process: the negative curvature that develops at the dent grows beyond the 
termination point of the upper branch in Fig. 2a and a transition to the lower branch takes 
place. This results in the reversal of propagation direction and the nucleation of a spiral- 
wave pair. The positive slopes of the upper and lower branches ai k — indicate transverse 
instabilities of the two planar Bloch front solutions. 

The proximity to the front bifurcation is essential for spontaneous spiral-wave nucleation. 
Farther from the bifurcation initial dents may grow due to the transverse instability but not 
nucleate spiral waves. This is demonstrated in Figs. 7 where the same initial conditions as 
in Fig. 6a are chosen. The initial almost planar front develops a dent but the dent retracts 
rather than nucleate a spiral-wave pair. 
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Fig. 6. A numerical solution of the kinematic equations (l)-(3) demonstrating the nucleation of a 
spiral- wave pair from an unstable propagating front. Left column: the Co(s) and k{s) profiles. Right 
column: the front Une shape in the x — y plane. A small perturbation on the initially near planar 
front (a) grows (b) and nucleates a pair of spiral waves (c). The front regions in the left frame of 
(d) each correspond to the core of the rotating spiral in the right frame. Parameters: ai = 4.0, 
ao = -0.0001, e = 0.0115, 6 = 1.063. (a)-(d) are at t = 0, 116, 136, 142. 

5 Conclusion 

We have developed a new set of kinematic equations for front dynamics in nearly symmetric 
bistable media. The equations are quantitatively valid for slowly propagating, weakly curved 
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Fig. 7. A numerical solution of the kinematic equations (l)-(3) with the same initial conditions as 
in Fig. 6 but farther away from the NIB bifurcation. Left column: the Co(s) and k{s) profiles. Right 
column: the front line shape in the x — y plane. The transverse instability causes a small dent (a) 

to grow (b) but since the system is far enough from the NIB bifurcation, no spiral-wave nucleation 
occurs (c), and the dent contracts (d). Parameters: ai = 4.0, oq = —0.0001, e = 0.0115, 6 = 1.055. 
(a)-(d) are at t = 0, 175, 220, 240. 

fronts. They describe the growth of transverse perturbations, the core structure of spiral 
waves, and the dynamics of spiral-wave nucleation. Away from the NIB bifurcation, where the 
front velocity is no longer a slow degree of freedom, the algebraic Cn-n relation is recovered. 
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The process of spiral-wave nucleation involves local transitions between the counterprop- 
agating Bloch fronts. In the present work these transitions were driven by curvature per- 
turbations. Front transitions and spiral nucleation can also be driven by other intrinsic 
perturbations, such as nonlocal front interactions or interactions with boundaries [16]. Such 
interactions, however, have not been included in the present derivation. They are excluded 
by the choice of the boundary conditions in Eqns. (17). Front interactions are not important 
for the study of symmetric spirals (Figs. 3) or the initial nucleation of a spiral-wave pair 
from a planar front (Figs. 6). They do, however, affect nonsymmetric spiral- waves (oq 7^ 0) 
and might play an important role in the meander instability of a spiral tip [41-43] . Front in- 
teractions are also essential for the formation of labyrinthine patterns in the Ising regime [7] . 

The kinematic equations generalize an earlier approach based on a geometric equation for 
curvature supplemented by a linear normal velocity - curvature relation [22-25]. In that 
case the speed of a planar front is taken to be a constant, determined by the parameters 
of the system, and no distinction is made between the two types of Bloch fronts. This 
approach has been applied to traveling waves in excitable media modeling a pulse stripe as a 
single curve [44-46]. It has also been applied to spiral- wave dynamics with phenomenological 
assumptions to describe the motion of the spiral tip, the free end of the curve [22-25,34]. No 
phenomenological assumptions are needed in applying the generalized kinematic equations 
as they naturally describe the core (or tip) structure of a spiral-wave. 

Finally, we note that the two-dimensional spiral-wave nucleation problem in the original 
bistable medium is reduced to a one-dimensional "droplet" nucleation problem in the kine- 
matic equations [47,48] . This may simplify the evaluation of the critical curvature perturba- 
tion required to nucleate spiral waves. 
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